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The  design  and  operation  of  Beam- 
con  III,  the  latest  linearity  measure- 
ment instrument  using  the  beam 
addition  method  in  the  detector  metrol- 
ogy program  at  the  National  Institute 
of  Standards  and  Technology,  is  de- 
scribed. The  primary  improvements  in 
this  instrument  are  the  reduction  of 
stray  radiation  to  extremely  low  levels 
by  using  three  well-baffled  chambers,  a 
larger  dynamic  range,  and  an  additional 
source  entrance  port.  A polynomial  re- 
sponse function  is  determined  from  the 


data  obtained  by  this  instrument  using 
a least-squares  method.  The  linearity  of 
a silicon  photodiode-amplifier  detector 
system  was  determined  to  be  within 
0.054  % (2cr  estimate)  over  nine 
decades  of  signal. 
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1.  Introduction 

The  underlying  assumption  in  the  use  of  most 
detection  methods  for  the  measurement  of  optical 
radiation  is  that  the  detector  output  signal  is  di- 
rectly proportional  to  the  input  radiation  flux.  This 
proportionality  is  defined  as  linearity  and  con- 
versely the  departure  from  proportionality  is  de- 
fined as  nonlinearity.  Optical  radiation  detectors 
are  typically  converters  of  optical  power  (or  photon 
flux)  to  measurable  electrical  parameters,  such  as 
current,  voltage,  or  pulse  frequency.  Therefore,  in 
radiometric  applications,  the  linearity  of  a detec- 
tor’s signal  is  not  only  a function  of  the  efficiency 
of  the  physical  radiation  detector  (transducer),  but 
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also  of  the  peripheral  electronics.  A well  known 
example  of  this  is  a photomultiplier/photon-count- 
ing system  where  significant  nonlinearity  is  induced 
at  high  count  rates  by  the  dead  time  of  the  pulse 
counter.  All  detector  systems,  consisting  of  detec- 
tors and  peripheral  electronics,  are  characterized 
by  some  degree  of  nonlinearity.  The  exact  positions 
of  these  nonlinear  regions  are  functions  of  both  the 
detector  and  the  accompanying  electronics  and 
should  be  determined  experimentally  for  measure- 
ments of  the  highest  accuracy.  There  are  essentially 
two  practical  choices  when  confronted  by  nonlinear 
behavior  of  a detector  system:  first,  to  view  it  as  a 
performance  limit,  a maximum  deviation  from  lin- 
earity over  a given  range  of  input;  second,  to 
correct  the  measured  signal  given  the  known 
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nonlinearity  of  the  detector.  This  latter  solution 
presupposes  that  the  detector’s  nonlinearity  is  well 
characterized  and  should  be  determined  with  the 
measurement  electronics  that  are  to  be  routinely 
used  with  the  detector. 

Concern  about  detector  linearity  and  the  devel- 
opment of  instruments  for  its  measurement  has  a 
long  history  at  the  National  Institute  of  Standards 
and  Technology,  formerly  the  National  Bureau  of 
Standards  (NBS).  Sanders  [1]  reported  in  1972  the 
results  of  photocell  linearity  measurements  using  a 
multi-aperture  design.  Also  in  1972,  Mielenz  and 
Eckerle  [2]  utilized  a double-aperture  method  to 
characterize  the  NBS  Reference  Spectrophotome- 
ter. In  1984  Saunders  and  Shumaker  [3]  reported 
measurements  using  an  automated  radiometric  lin- 
earity tester.  This  instrument  was  called  a beam- 
conjoiner. 

Based  on  this  initial  work,  we  now  report  on  our 
latest  generation  of  beamconjoiner  technology, 
Beamcon  III.  This  instrument  was  specifically  de- 
signed to  characterize  radiometric  detectors,  par- 
ticularly at  low  flux  levels  approaching  the  noise 
equivalent  power  (NEP)  of  the  detector.  To  facili- 
tate low  level  measurements,  this  instrument  has 
three  well-baffled  chambers  to  reduce  the  stray  ra- 
diation to  extremely  low  levels.  Other  improve- 
ments are  additional  filters  to  give  a larger  dynamic 
range  and  a second  source  entrance  port  for  a laser 
source  or  for  source  mixing.  The  data  analysis  al- 
gorithm of  Beamcon  III  calculates  the  responsivity 
of  the  detector  system  assuming  a polynomial  rela- 
tionship between  the  input  flux  and  the  output  sig- 
nal from  the  detector  system.  Beamcon  III  will 
serve  as  an  automated  instrument  for  routine  de- 
tector linearity  characterization  and  for  determin- 
ing the  linearity  of  NIST’s  monochromator-based 
systems  (e.g.,  spectroradiometers  and  spectropho- 
tometers). 

2.  Description  of  Beamcon  III  and  its 
Operation 

A schematic  diagram  illustrating  the  technique 
of  beamconjoining  is  shown  in  Fig.  1.  The  input 
beam  from  the  source  is  split  into  two  optical  paths 
by  the  first  beamsplitter  BSl  and  variably  attenu- 
ated by  filters  on  either  wheel  W1  or  W3.  After 
reflection  off  mirrors  Ml  and  M2,  the  two  beams 
are  combined  by  beamsplitting  BS2  and  attenuated 
again  by  a third  filter  on  wheel  W2  before  falling 
on  the  detector.  All  filters  are  metallized  neutral 
density  filters  with  quartz  substrates,  and  the  filter 


W3 


Fig.  1.  Schematic  optical  diagram  of  Beamcon  III  showing  the 
primary  optica!  elements.  BS;  beam  splitter;  W:  filter  wheel; 
and  M:  mirror. 

wheels  are  tilted  to  preclude  interreflection  be- 
tween filters.  In  addition  the  beamsplitters  are 
wedged  to  prevent  etaloning. 

A complete  schematic  diagram  of  Beamcon  III, 
including  all  ports,  optics,  and  chambers,  is  shown 
in  Fig.  2.  Beamcon  III  provides  four  basic  advan- 
tages over  the  previous  technology:  1)  each  cham- 
ber is  optically  isolated  from  the  next  chamber 
except  for  a 3 inch  diameter  throughput  port, 
thereby  reducing  scattered  radiation;  2)  there  are 
two  source  input  ports  for  both  laser  and  broad- 
band sources  either  singly  or  concurrently  for 


Chamber  3 


Fig.  2.  System  layout  of  Beamcon  III.  Chamber  1 collimates  the 
beam  from  the  source.  Chamber  2 splits  the  beam  into  two 
paths  and  recombines  it.  Chamber  3 has  a spatial  filter  SF  to 
eliminate  stray  light  and  focuses  the  beam  onto  the  detector. 
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source  mixing;  3)  a steel  base  plate  allows  for  mag- 
netic mounting  of  additional  optical  components; 
4)  an  external  translation  stage  on  the  broadband 
input  port  can  vary  the  numerical  aperture  and  col- 
limate the  source  beam.  A microcomputer  controls 
the  positions  of  the  filter  wheels  and  records  the 
signal  from  the  detector. 

An  incoherent  light  source  (Source,  Fig.  2),  ei- 
ther a quartz  halogen  or  arc  lamp,  is  placed  on  a 
translation  stage  in  front  of  the  input  port  to 
Chamber  1.  The  input  power  of  the  source  beam  is 
varied  either  by  changing  the  aperture  of  the  iris 
diaphragm  (Iris  1)  or  by  translating  the  Source,  or 
both.  The  input  beam  in  Chamber  1 is  folded  by  a 
flat  mirror  (Mf)  and  collimated  by  a concave  spher- 
ical mirror  (Me),  which  projects  the  beam  through 
the  input  port  of  Chamber  2.  For  narrowband  op- 
eration with  an  incoherent  light  source,  an  interfer- 
ence filter  can  be  placed  at  the  input  port  of 
Chamber  2.  The  first  beamsplitter  (BSl)  yields  two 
separate  beams  which  are  independently  attenu- 
ated by  two  filter  wheels  (Wl)  or  (W3).  These  fil- 
ter wheels  contain  four  neutral  density  filters  with 
differing  optical  densities.  The  two  beams  are 
folded  by  mirrors  (Ml)  or  (M2)  and  recombined 
into  a single  beam  by  beamsplitter  (BS2).  This 
beam  is  further  attenuated  by  passing  through  a 
third  filter  wheel  (W2)  containing  five  different 
neutral  density  filters.  The  beam  then  passes  into 
Chamber  3,  which  contains  a set  of  concave  mirrors 
(M3  and  M4).  A spatial  filter  (SF)  at  the  focal 
point  between  the  two  mirrors  filters  out  scattered 
and  diffracted  light,  resulting  in  a uniform  output 
field  at  the  detector.  A concave  mirror  (M5)  fo- 
cuses the  beam  onto  the  detector.  Translational 
shifting  of  M5  controls  the  image  plane  of  the  in- 
strument and  varies  the  irradiance  at  the  detector. 
The  detector  is  generally  overfilled  to  prevent  arti- 
facts due  to  different  cut-off  edges  in  the  optical 
path  (optical  stop  variation)  and  to  the  local 
nonuniformity  of  the  detector.  In  addition  environ- 
mental parameters  such  as  humidity  and  tempera- 
ture are  constantly  monitored  at  the  detector.  A 
monochromator  can  also  be  attached  at  the  exit 
port  to  select  a wavelength  for  narrowband  opera- 
tion. 

Beamcon  III  has  a second  input  port  into  Cham- 
ber 2,  which  can  be  used  singly  or  in  concert  with 
the  beam  from  Chamber  1.  Generally  the  input 
port  of  Chamber  2 is  used  for  sources  such  as  He- 
Ne  or  Ar"^  lasers.  A beam  expander  BE  is  used  to 
minimize  scattering  losses  and  filter  transmittance 
nonuniformity.  The  laser  can  also  be  used  to  align 
both  optical  paths  by  removing  the  reflecting  mir- 


ror (Mr)  and  inserting  the  mirror  (Ma)  to  guide  the 
laser  beam  through  the  center  of  the  lamp  filament 
and  of  every  optical  element.  The  mixed  source 
method  has  not  been  attempted  as  yet,  but  is 
planned  in  future  studies,  as  it  has  advantages  for 
studying  differential  spectral  detector  nonlineari- 
ties [4]. 

Stepping  motors  control  the  filter  wheel  posi- 
tions, each  with  400  steps/revolution.  A tolerable 
error  is  one  step,  therefore  the  error  in  filter  posi- 
tion is  0.9°.  The  space  between  two  filters  on  each 
wheel  acts  as  a shutter  position  for  a dark  signal 
measurement.  Thus,  the  total  number  of  filter  com- 
binations is  150,  and  a full  150  measurements  is 
called  a pass.  Of  these,  30  are  dark  signal  measure- 
ments: 5 when  both  Wl  and  W3  are  shuttered  and 
25  when  W2  is  shuttered.  The  remaining  120  filter 
combinations,  yielding  120  response  signal  mea- 
surements, are  randomized  in  each  pass  to  prevent 
hysteresis.  The  dark  signal/filter  combinations  are 
randomized  and  inserted  after  every  five  measure- 
ments. Each  dark  signal  is  used  to  correct  the  suc- 
ceeding five  response  signal  measurements. 
Generally,  five  passes  are  taken  and  averaged  to 
reduce  the  random  error  and  to  collect  information 
about  possible  source,  filter,  and  detector  drift. 
The  maximum  throughput  is  7 % at  a wavelength 
of  632.8  nm,  and  the  total  attenuation  factor  is 
about  3000  against  a light  source. 

The  output  beam  is  partially  vertically  jX)larized 
for  an  unpolarized  source  since  each  tilted  surface 
preferentially  reflects  the  vertical  component  of 
the  incident  beam.  The  ratio  of  the  polarization  in 
the  vertical  direction  to  that  in  the  horizontal  di- 
rection is  1.17.  Thus,  it  is  critical  to  keep  the  polar- 
ization direction  of  a source  (e.g.,  laser)  constant  in 
order  to  eliminate  this  polarization  preference. 

3.  Polynomial  Response  Function  Fitting 

The  principle  of  operation  of  Beamcon  III  is  the 
beam  addition  method,  in  which  different  filter 
wheel  combinations  result  in  different  fluxes  at  the 
detector.  This  method  assumes  that  the  flux  from 
the  source  remains  constant  during  all  measure- 
ments and  that  the  output  flux  at  the  detector  is 
the  sum  of  the  two  fluxes  along  the  two  different 
paths. 

The  final  flux  (f)  at  the  detector  is  the  sum  of  the 
fluxes  (f)i  and  (f>2  from  the  two  paths.  Indexing  the 
fluxes  by  the  filter  combination. 
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where  i equals  0 to  4 designates  the  filter  on  wheel 
Wl,y  equals  0 to  4 the  filter  on  W3,  and  k equals  0 
to  5 the  filter  on  W2.  In  each  case,  an  index  value 
of  0 corresponds  to  an  opaque,  shuttered  position. 
For  a single  path,  there  are  30  possible  filter  com- 
binations. Of  these,  no  flux  is  transmitted  for  10 
combinations.  Thus,  each  path  has  20  filter  combi- 
nations that  result  in  a response  signal. 

Since  the  detector  is  assumed  to  be  nearly  linear, 
it  is  appropriate  and  convenient  to  assume  an  n -th 
order  polynomial  response  function  for  the  detec- 
tor [3],  so  that 

(l){iJ,k)=ro  + ris(ijJc)+ris\ijJc)  + ...  + r„s''{ijJc), 

(2) 

where  s(ijjc)  is  the  measured  response  signal,  cor- 
rected by  the  dark  signal  and  indexed  by  the  filter 
combination,  and  ro,  n,  ri,  ...,  r„  are  the  response 
function  coefficients.  The  coefficient  ri  is  set  equal 
to  one,  both  because  the  detector  is  assumed  to  be 
nearly  linear  and  to  normalize  the  flux  to  the  same 
units  as  the  signal. 

Combining  Eqs.  (1)  and  (2),  with  n = 1,  yields 

(l)i(iJc)  + (fy2(jJc)=ro  + s(ijJc)+r2S  \i  J Jc  ) 

+ ...  + rnS''{iJJc).  (3) 

There  are  120  expressions  of  the  form  of  Eq.  (3) 
from  the  120  distinct  response  signals  measured  in 
each  pass,  with  40 + n unknown  variables  from  the 
20  filter  combinations  along  each  path  and  the  n 
response  function  coefficients.  The  fluxes  of  shut- 
tered filter  combinations  are  set  equal  to  zero. 
These  120  expressions  are  solved  using  the  linear 
least-squares  technique  to  determine  the  values  of 
the  unknowns.  While  the  values  of  the  polynomial 
response  function  coefficients  are  the  primary  goal 
of  the  measurements,  the  values  of  the  fluxes  are 
useful  for  monitoring  the  stability  of  the  filters. 

Solving  for  both  the  fluxes  along  the  paths  and 
the  response  function  coefficients  provides  two 
ways  in  which  to  determine  the  total  flux  at  the 
detector,  using  Eqs.  (1)  and  (2),  and  therefore 
serves  to  verify  the  results.  A measurement  run 
yields  120  residuals  between  the  total  fluxes  calcu- 
lated in  these  two  ways.  If  there  is  a slope  in  a 
time-ordered  plot  of  these  residuals,  then  either 
dark  signal  or  source  intensity  variations  likely  ex- 
ist. Also,  the  standard  deviation  of  these  residuals 
indicates  the  noise  level  present  during  the  mea- 
surement, either  from  the  source,  filters,  or  detec- 
tor. For  most  detectors,  a second-order  polynomial 


accurately  fits  the  data.  However,  if  the  detector  is 
highly  nonlinear,  the  power  of  the  polynomial  is 
increased  until  the  residuals  are  normally  dis- 
tributed about  zero.  If  the  response  function  coeffi- 
cients remain  unreasonable,  then  either  a 
polynomial  is  not  an  appropriate  fitting  function  or 
the  detector  and  its  associated  electronics  are  not 
functioning  properly. 

The  algorithm  for  determining  the  40 -I- n un- 
known variables  from  a complete  measurement  was 
extensively  tested  using  synthetic  data  sets  with 
known  response  function  coefficients  for  second- 
order  polynomials.  It  was  also  tested  using  syn- 
thetic data  sets  corresponding  to  two  cases  of 
nonlinear  behavior,  supraresponsivity  and  satura- 
tion at  high  flux  levels  [5].  The  algorithm  was  suc- 
cessful in  fitting  these  two  behaviors  with  a third- 
and  a sbeth-order  polynomial,  respectively. 

4.  Linearity  Measurement  of  a Silicon 
Photodiode 

Beamcon  III  was  used  to  determine  the  linearity 
of  a silicon  photodiode-amplifier  detector  system 
designed  and  constructed  at  NIST  for  multi-decade 
performance  and  described  previously  [6].  The  lin- 
earity was  determined  for  amplifier  gains  from  10“ 
to  10'°.  The  broadband  source  was  a 1000  W tung- 
sten-halogen lamp  operated  at  a constant  current 
of  7.6  A.  In  order  to  cover  the  entire  signal  range 
of  the  detector  system  without  changing  the  spec- 
tral distribution  of  the  source,  the  flux  from  the 
lamp  was  varied  either  by  changing  the  numerical 
aperture  of  Iris  1 in  Fig.  2 or  by  placing  neutral 
density  filters  in  Chamber  1,  or  both. 

The  relative  responsivity  of  a detector  system  in- 
dicates deviations  from  linearity.  Once  the  re- 
sponse function  coefficients  have  been  determined, 
the  total  normalized  flux  for  a given  signal  is  given 
by  the  response  function  [the  right-hand-side  of 
Eq.  (3)].  If  the  detector  system  is  linear,  the  total 
normalized  flux  is  simply  the  signal.  Therefore,  the 
relative  responsivity  at  a given  signal  is  the  re- 
sponse function  divided  by  the  signal.  The  results 
from  one  measurement  of  the  linearity  of  the  de- 
tector system  at  a gain  of  10*  is  shown  is  Fig.  3, 
where  the  relative  responsivity  is  plotted  as  a 
function  of  signal.  The  response  function  was  a 
second-order  polynomial  with  coefficients 
/-o= -5.53  X 10“’  and  /■2=  - 1-56  x 10“*.  Over  the 
three  decades  of  signal  measured  at  this  gain,  the 
detector  system  is  linear  to  within  0.027  % (2cr 
estimate). 
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Fig.  3.  Relative  responsivity  as  a function  of  signal  for  a silicon 
diode-amplifier  detector  system  at  a gain  of  10^.  The  error  bar 
indicates  the  Zo- estimated  uncertainty  of  the  relative  responsiv- 
ity. 

The  maximum  signal  for  the  detector  system,  de- 
termined by  the  source  and  throughput  of  Beam- 
con  III,  was  about  10“^  A.  At  a bandwidth  of  0.67 
Hz,  the  NEP  of  the  detector  system  resulted  in  a 
minimum  measurable  signal  of  about  10“’^  A.  In 
order  to  characterize  the  detector  system  between 
these  limits,  linearity  measurements  were  per- 
formed at  all  seven  amplifier  gain  ranges.  A sec- 
ond-order polynomial  response  function  was 
calculated  for  each  range,  and  the  relative  respon- 
sivity over  nine  decades  of  signal  is  shown  in  Fig.  4. 
The  detector  system  is  linear  to  within  0.054  % (la- 
estimate)  over  the  entire  range,  and  to  within 
0.209  % (2o- estimate)  for  currents  less  than  10"" 
A.  Thus,  Beamcon  III  was  able  to  determine  the 
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Fig.  4.  Relative  responsivity  as  a function  of  signal  for  a silicon 
diode-amplifier  detector  system  over  its  entire  signal  range  of 
nine  decades.  The  error  bars  indicate  the  estimated  2<r  uncer- 
tainty of  the  relative  responsivity  for  signals  less  than  10"  A (on 
the  left)  and  over  the  entire  range  of  signal  (on  the  right). 


linearity  of  this  detector  system  for  signals  ap- 
proaching its  NEP. 


5.  Conclusions 

The  design  and  operation  of  our  latest  genera- 
tion of  beamconjoining  technology,  Beamcon  III, 
was  detailed.  Using  this  instrument,  the  linearity  of 
a silicon  photodiode-amplifier  detector  system  was 
determined  over  nine  decades  of  signal,  approach- 
ing its  NEP  at  the  smallest  signals.  Beamcon  III 
will  permit  both  routine  calibration  of  detector  sys- 
tems and  research  into  detector-amplifier  effects 
on  nonlinearity.  Future  studies  on  the  linearity  of 
solid  state  detector-amplifier  systems,  photomulti- 
plier photon-counting  systems,  and  multichannel 
detectors  (i.e.,  diode  arrays  and  charge  coupled 
devices)  are  planned,  as  well  as  on  the  capabilities 
of  Beamcon  III  for  source  mbcing. 
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